Neutron localization induced by the pairing field in an inhomogeneous neutron matter 
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It is shown that in an inhomogeneous neutron matter the pairing field bounds neutrons around 
the Fermi level leading to the formation of localized Andreev states. In the case of the inner crust 
of neutron stars the localization length has been determined as a function of the nuclear density. 
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The inhomogeneous nuclear matter forms the inner 
crust of neutron stars at densities ranging from lO""* 
fm~'^ to 10~^ fm~'^. The nuclear matter is thus diluted 
and the nucleons are strongly paired through the attrac- 
tion in the s-wave channel. Since the proton fraction does 
not exceed a few percent in this density region the main 
contribution to the pairing energy comes from neutrons. 
The dependence of the pairing field on the neutron den- 
sity has been a subject of an intensive theoretical studies 
(see [H and references therein). Approaches based on 
the BCS theory predict the maximum of the s-wave pair- 
ing gap A to be around SMeVQ, whereas the inclusion 
of polarization corrections decreases this value to about 
IMeyQ. Recently, the non-perturbative calculations for 
fermions in the unitary limit predict the pairing gap at 
T = of about 0.6eF, where is the Fermi energy Q. 
Since the neutron matter at the densities between 10"'^ 
fm~'^ and 10~^ fm~^ is close to the unitary limit, it sug- 
gests that the pairing gap is a considerable fraction of 
the Fermi energy. 

The proton admixture drives the nuclear matter to- 
wards an inhomogeneous state and results in the appear- 
ance of nuclear impurities immersed in the superfluid 
neutron liquid. Inside impurities the density reaches the 
nuclear saturation density. The lattice of nuclei is sta- 
bilized by the Coulomb interaction and becomes denser 
as the average density increases, and eventually melts 
due to the interplay between Coulomb, surface and shell 
energies @, 0] ■ Consequently, the inhomogeneous neu- 
tron density distribution induces the spatial variations 
of the pairing potentialjl, Q. Because the structure of 
the major part of the crust is predicted to be periodic, 
the pairing field is periodic as well. Moreover pairing 
at nuclear saturation density is usually weaker than in 
the dilute neutron matter (except for a very dilute sys- 
tem with a density smaller than 10~'^fm~'^). Hence the 
local minima of the pairing field coincide with the posi- 
tion of nuclear impurities. The scale of spatial variations 

of the pairing field is governed by the coherence length 
1 ep 

^ — — , where kp stands for Fermi momentum. Al- 

though the pairing gap is large, it is still smaller than ep 
and thus the variations of the pairing field are smoother 
than that of the nuclear mean-field potential. 



The inhomogeneity of the pairing field induces an ad- 
ditional scattering of quasiparticles. Because the pairing 
field is periodic, an additional Bragg scattering gives rise 
to corrections to the neutron band structure. Therefore 
the natural question appears to be whether this effect is 
large, or can be neglected. This article provides an an- 
swer to this question and discuss the impact of the inho- 
mogeneity of the pairing field on the quasiparticle states. 
It is organized as follows: first, I describe the mecha- 
nism of bound state formation in the schematic model. 
Second, it will be shown, that contrary to the normal 
three-dimensional potential, in the case of a pairing well, 
there is always a bound state, no matter how shallow 
the well is. Finally, I will present the variational calcu- 
lations for realistic pairing potentials through the inner 
crust and estimate the density region where the localized 
states appear. 

In order to investigate the problem one must realize 
that the scattering of quasiparticles on the pairing field 
is qualitatively different than that on the normal poten- 
tial. Namely, the scattering matrix contains in this case 
the normal part which comprises the scattering ampli- 
tude of hole into hole and particle into particle, but in 
addition it contains the amplitude of scattering of par- 
ticle into hole and vice versa. This, so called Andreev 
reflection has an unusual property that the particle/hole 
is reflected almost exactly in the same direction as the 
incoming hole/particle [lO|. In fact, it can be shown that 
the Andreev reflection dominates around the Fermi level 
and gives rise to the formation of the so-called Andreev 
states [13, in I- 

The spatial variations of the pairing field may be sim- 
ulated in the first approximation by the simple model 
consisting of the spherical pairing well: 
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and|Ao„t| > |Ai„|. In the absence of superflow the phase 
of the pairing field is constant and can be chosen to be 
zero. The BdG equation takes the form: 
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where T is the kinetic energy operator (for simplicity we 
disregard the normal potential) , E is the quasiparticle en- 
ergy, u{r) and v(r) are the particle and hole amplitudes of 
the quasiparticle wave function. In the simplest Andreev 
approximation, which will admit the analytic solution, 
the wave function is decomposed into the rapidly oscil- 
lating part, and the smooth part, related to the changes 
of the pairing field [13] : 
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= cxp(ikF • r) 
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where u(r) and v{r) changes slowly within the scale deter- 
mined by the coherence length. Neglecting terms of the 
order of l/{kp^) the original equations transform into: 
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We have chosen units: ^ — 
kp 



v{r) J y v{r) 

(4) 

1, i.e. T = -V^ and 
fi. Although the Andreev approximation must be 
used with care for a very dilute neutron matter, it is in- 
structive to study this limit (i.e. kp^ — > c». A/ ^ 0), as 
the consequences of this oversimplified approach remain 
valid for a more realistic cases. 

The astonishing implication of eq. (|4]) is that the three- 
dimensional problem is transformed into a one dimen- 
sional, and for each path specified by the vector kp ex- 
ists a solution. This result is a direct consequence of the 
focusing property of the Andreev reflection [lol 11 1. In 
order to investigate the L = state, one must consider 
the quantization condition along the diameter: 



A(0+,(/)_)exp(2igi?) = 1, 
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The closer inspection of the above result reveals that the 
equation ([S]) always admits a solution for E in the interval 
(Ai„, Ao„(), no matter how small is the difference |Ai„| — 
\/S.out\- Namely, in the hmit |Ai„| < |AoMt| — > the 
quantization condition becomes: 



exp(2i(gi? + $)) = 1, 



(7) 



where 



cos $ = 



y = 
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Since qR — > and y < x < 1 it is easy to find that the 
solution reads: 
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Another limit relevant for the neutron star crust is when 
A,;„ 0. In that case one gets: 

1 E E 

— kpR arccos -7 — Tim, m = 0, ±1, ±2, ... (10) 

2 Aoiif 

The ^xamination of both limits reveals that once 

kpR — 2^ and are fixed, then — is constant as 

M Aout Aotit 

well. This property is worth mentioning in the context of 

the neutron star crust, because there i? oc ^, as the size of 

nuclear impurities is small with respect to the coherence 



length. In this case clearly kpR- 
E 



A„ 



is approximately 



constant and thus 



stays fixed as well, through the 



whole inner crust. Consequently, one may expect that 
the size of the tail of the wave function bound by the 
pairing field, 7], will be proportional to ^ and r//^ will be 
approximately constant through the crust. In the follow- 
ing we will show that this is indeed the case. 

Let us turn now to the three-dimensional description, 
in order to study whether the above results survive in a 
more realistic case. In the variational formulation of the 
problem ([2]), one must consider the expectation value of 
the square of the energy: 



{E') = / d^rx+(r)r!+r!x(r), 



where 
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and T = — and k\ = ji. We use the ansatz so that 
the wave function x(r) can be decomposed similarly to 
©: 



X(r) = 0fc^(r)^^fcjr), 



(13) 



where tj^k (r) is a spinor, and the following conditions 
hold: 



V20fc^(r) = 40fcF(i-), 
-VVfc,(r) « fc'V'fcJr), 



(14) 
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Because the problem is time-reversal invariant then x(r) 
can be chosen to be real, and consequently the term 
([V^, A(r)]} ~ 0. Hence ([TT|) can be rewritten as: 



where the first term can be estimated as follows: 



(15) 
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with a, (3 being numerical factors of the order of unity. 
In order to estimate {V"^) let us assume now that A(r) = 
Ao„t + 5A(r) and (5A(r) ^ for r < i?. If wc treat 
(5A(r) as a small perturbation on top of the homogeneous 
pairing field Ao„t then: 



^out + 2Ao„f 
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where JA — 5A(ro) and ro < R. Clearly, the last in- 
tegral contains a smooth and a rapidly oscillating part. 
The smooth part can be extracted using the mean value 
theorem: 



r<R 



d'rx' (r)x(r) = ^k, (ri)^fe, (ri) / d^^k, (r)) 
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where ri < R. If kpR ^ 1 then the right hand side 
integral is proportional to R instead of R'^ as one would 
expect. This fact can be understood in the classical pic- 
ture since the probability of finding a free particle inside 
some area is proportional to its linear size. On the other 
hand, since the above integral is dimensionless it should 
be a function of k,^R. Hence one may write: 



2AoutSAf{k^R), 



(19) 



where /(k^iR) — s- jk^/R, if A:^ — s- and 7 > is a nu- 
merical constant. Using the above estimates one may 
express (E"^) as a function of k^. Because we would like 
to answer the question whether the minimum of energy 
corresponds to fc^ > 0, it is sufficient if we consider the 
limit of fc,, being very small. In this case: 



7 2 2 
^k^^a 



Al^t + 2Aout5A-fk^R. (20) 



Note the striking similarity with the estimation for the 
ground state energy of a normal potential in one dimen- 
sion. The condition for the minimum of {E"^) gives the 
relation for /c„: 
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(21) 



pairing potential well, whereas a positive value would be 
equivalent to a potential barrier. 

Hence the scattering on the pairing potential has a 
more severe effect on the quasiparticle spectrum as the 
pairing field turns effectively the three-dimensional prob- 
lem into a one-dimensional, and thus the scattering is 
more effective. 
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FIG. 1: The results of minimization of (-B^) with respect 
to localization length 77 (see text) plotted as a function of the 
nuclear density in the inner crust of neutron star. The dashed 
line corresponds to the BCS pairing and the solid line denotes 
the pairing with inclusion of polarization corrections. The 
lower subfigure shows the coresponding pairing gaps, used in 
the calculations. 

In order to quantitatively study this effect we minimize 
(Jii^) in a class of wave functions of the form: 

expf-fc^r^l { " 

V 



(22) 



with/dVx-(r)x(r) = 1. The pairing field was assumed 
to have the form: 

A(r) = A,,,t + (5Aexp(-(r/0'), (23) 



It shows that (E^) is minimized by kr, > 0, if A^^t^A < 0. '^^'^^'^ " ^"^ " ™^ ^ 

The last condition is obvious as it corresponds to the 



In 

- — . The actual 

kp Aout 

parameters Aout , Ai„ requires some knowledge about the 
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density dependence of the pairing field. The two sets 
of calculation has been performed. In the first one, 
the pairing field ^out was assumed to have the den- 
sity dependence obtained from the BCS approach with 
Argonne vl8 interaction[2|, whereas Ai„ corresponds 
to the neutron pairing at the nuclear saturation den- 
sity. The second set of calculation was performed for 
pairing extracted from the s-wave N-N scattering data 

A « - — — ^ exp — — - — - , which is sup- 

\e) 2m V 2tan(5(fcF)/ 
posed to include the polarization corrections (JJIJ. The 
results are presented in the Fig. 1. In the upper sub- 
figure the Rc/v is plotted, where rj — and Rc is 
the Wigner-Seitz cell taken according to the calculations 
from refs. [l^ . 

The results indicate a strong dependence of rj on the 
pairing field. In the case of stronger pairing the quasipar- 
ticle states become practically localized for lower densi- 
ties, i.e. Rc/fJ > 1. For weaker pairing the impact of the 
inhomogeneous pairing field is less pronounced although 
the ratio Rc/rj reaches 0.6. It has to be noted that the 
ratio ^/{£P)/Aout ~ 0.9 throughout the crust. Conse- 
quently 77/^ « 50. It supports the conlcusion based on 
the schematic model. 

Note that in the calculations the mean nuclear poten- 
tial was set to zero. Otherwise the results would have 
been obscured by the existence of an additional scatter- 
ing which would alter the rapidly oscillating part of the 
wave function. The current approach corresponds to the 
Thomas-Fermi approximation for the mean-field poten- 
tial. Thus the results represent the average behavior of 
the quasiparticle states around the Fermi level. The so- 
lution of the full BdG or HFB equations can change the 
above results in the following way: around the Fermi 
level will appear resonances induced by the mean-field. 
However their pairing properties will be practically unal- 
tered, as compared to our model, unless their energy will 
deviate from the Fermi level more than |Aout|. 

The consequences of the presented results are twofold. 
First, they suggest that the scattering of quasiparticles 
has to be taken into account if the neutron band struc- 
ture is to be determined. Consequently, for this purpose 
the BCS approach is too poor, as there the quasiparticle 
scattering is neglected. On the other hand BCS is suf- 
ficient if one is interested in the bulk properties of the 
inner crust: energy density, density distribution, proton 
to neutron fraction, etc. It is interesting to note that the 
need for HFB (or BdG) approach has been noticed in the 
context of nuclei close to the neutron drip line, long time 
ago . There the argument originated from the nuclear 
density localization requirement. Here, the main argu- 
ment comes from the realization that the pairing field 
exhibits spatial variations in the inhomogeneous nuclear 
matter. When the matter becomes diluted the scatter- 
ing on the pairing potential becomes more effective and 



contribute to the band structure around the Fermi level. 

Second, the consequence of neutron localization by the 
pairing field may turn out to be important for transport 
properties. Namely, in this case the effective mass has 
to be determined and thus the band structure has to be 
calculated Within the approach we presented here, 
one may estimate the dispersion relation for the band cor- 
responding to bound neutrons. Taking into account only 
the nearest neighbour interaction of the localized neu- 
trons one obtains: E{k) oc const. + exp{— ^{Rc /r])^)k^ ■ 
Note that in the limit of Rc/v > 1 the coefficient in 
front is practically zero and the effective mass tends to 
infinity. This corresponds to complete localization of neu- 
tron states at the Fermi level around nuclear impurities. 
Clearly the detailed calculations need to be performed in 
order to investigate this problem. 

Summarizing, it was shown that the inhomogeneity of 
the pairing field in the neutron matter leads to the for- 
mation of neutron bound states around the Fermi level. 
They appear at an arbitrarily weak deviation from the 
homogeneous pairing field and are the consequence of 
the focusing property of Andreev reflection. 
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